Spin Dynamics and Multiple Reflections in Ferromagnetic Film in Contact with 

Normal Metal Layers 
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Spin dynamics of a metallic ferromagnetic film imbedded between normal metal layers is studied 
using the spin-pumping theory of Tserkovnyak et al. [Phys. Rev. Lett. 88, 117601 (2002)]. The 
scattering matrix for this structure is obtained using a spin-dependent potential with quantum 
well in the ferromagnetic region. Owing to multiple reflections in the well, the excess Gilbert 
damping and the gyromagnetic ratio exhibit quantum oscillations as a function of the thickness of 
the ferromagnetic film. The wavelength of the oscillations is given by the depth of the quantum 
well. For iron film imbedded between gold layers, the amplitude of the oscillations of the Gilbert 
damping is in an order of magnitude agreement with the damping observed by Urban et al. [Phys. 
Rev. Lett. 87, 217204 (2001)]. The results are compared with the linear response theory of Mills 
[Phys. Rev. B 68, 0144419 (2003)]. 
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PACS numbers: 75.40.Gb 75.70.-i 76.60.Es 73.40.- 



I. INTRODUCTION 

Spin dynamics in multilayers consisting of a thin ferro- 
magnetic film (F) imbedded between normal metal films 
(N) has received much attention recently. Ferromag- 
netic resonance (FMR) experiments on N/F/N struc- 
tures show that imbedding produces an additional spin 
damping 1,2 . A possible explanation of these observations 
is that the localized spins in the fcrromagnet interact, via 
s-d exchange, with the conduction electrons in the nor- 
mal metal 3 . Due to this interaction, a transfer of angu- 
lar momentum takes place from the precessing ferromag- 
net to the conduction electrons in the adjacent N layers. 
The reactive torque generated by this transfer acts back 
on the ferromagnetic magnetization causing an increased 
linewidth and a shift of the ferromagnetic resonance ab- 
sorption curve. 

An interesting formulation of this mechanism has been 
given in recent work by Tserkovnyak et al 4 . Extending 
the scattering approach of parametric charge pumping by 
Brouwer 5 to spin pumping, these authors calculate the 
spin current pumped through the N/F and F/N contacts 
of the N/F/N structure. The modifications of the spin 
dynamics, due to the induced spin current, are expressed 
using the Landau-Lifshitz-Gilbert (LLG) equation 6,7 



dm IT dm 

— = - 7 mxH eff + amx — (1) 

where H e ff is the effective magnetic field, m is the 
unit vector in the direction of the ferromagnetic magne- 
tization, 7 is the gyromagnetic ratio, and a is the Gilbert 
damping constant. 

According to Ref. [4], the spin pumping produces 
renormalization of the parameters 7 and a given by 
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where parameters cyq and 70 are the bulk values of a 
and 7. 5l is the Lande' factor, is the Bohr mag- 
neton, and M is the total moment of the ferromagnetic 
film. The superscripts (L) and (R) correspond to the left 
(N/F) and right (F/N) interfaces of the N/F/N structure. 

The parameters A\ L ' R ' and Ar L,R ^ stem from the spin 
currents pumped into the (L,R) normal layers (leads). 
These parameters are expressed in terms of the scatter- 
ing amplitudes as follows 4 
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) and {t^ n , t^ n ) are the reflection and 



transmission amplitudes for electrons with up and down 
spins. The expressions (4) and (5) are to be evaluated 
with transverse modes (m, n) taken at the Fermi energy. 

From Eqs. (2) and (3) we see, in conjunction with Eq. 
(1), that the parameter Aj is responsible for the shift 
of the FMR frequency, whereas A r increases the FMR 
linewith. For a symmetric N/F/N structure, we have 
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When the electrons scattered from one interface inter- 
fere incoherently with the other interface, the parameters 

A[ L ^ and A^ contribute independently of the parame- 
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ters A ( » R) and A r R) . 



As pointed out in Ref. 4, this applies 
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for ferromagnetic films that arc thick compared to the co- 
herence length A c ~ (k F — kp)^ 1 where k F and k F are the 
spin-up and spin-down Fermi wave vectors, respectively. 
In this case, the parameters Ai and A r can be expressed 
in terms of the mixing conductances of the F/N and N/F 
contacts 8 . 

This approximation suppresses the effect of multiple 
reflections in ferromagnetic layers that are not much 
thicker than A c . The physics is very similar to that of 
the Fabry-Perot interferometer in optics. One manifes- 
tation of such interferences is oscillation of interlayer ex- 
change coupling in F/N/F multilayers as a function of 
ferromagnetic-layer thickness. This effect has been first 
predicted from numerical calculations by Barnas 9 . Sub- 
sequently, Bruno 10 offered an explanation based on quan- 
tum interference picture. Results of his work have been 
confirmed by Bloemen et al. 11 who succeeded to see oscil- 
lations of the interlayer exchange coupling in Co/Cu/Co 
as the thickness of the Co layers is varied. 

The subject of the present work is to apply the spin- 
pumping approach to a model which allows multiple re- 
flections in the ferromagnetic layer to be manifested in 
the renormalization of the parameters a and 7. Primary 
motivation for this investigation comes from the recent 
work by Mills 12 . 

The starting point of Mills work is an approach 13 
which generalizes the Rudcrman-Kittel-Kasuya-Yosida 
(RKKY) theory to time-dependent source represented by 
the precessing vector m(t). In this approach, the renor- 
malization of the parameters 7 and a is found propor- 
tional to the frequency derivative of the real and imagi- 
nary parts of the transverse susceptibility, respectively. 

Mills calculates the transverse susceptibility from a 
model which allows for multiple reflections in the ferro- 
magnetic layer. Instead of representing the ferromagnetic 
layer by two independently scattering outside interfaces, 
(as done in Refs. [4] and [13 ]), the electrons in Ref. [12] 
are assumed to move in a spin-dependent potential of 
the trilayer in the form of a quantum well. This analysis 
yields an excess damping parameter a' = a — ao which 
exhibits strong quantum oscillations with the thickness 
of the ferromagnetic film. 

Ref. [13] and recent work of the present author 14 show 
that, for a ferromagnetic monolayer, the dynamic RKKY 
approach and the spin pumping theory yield the same 
results for the parameter a'. Thus, we anticipate that 
application of the spin-pumping method to scattering by 
a quantum well potential will produce renormalizations 
of the parameters a and 7 which also exhibit oscillations 
with the thickness of the ferromagnetic film. 

In the present investigation we confirm this anticipa- 
tion. However, it turns out that our results are, in several 
aspects, different from those of Ref. 12. Primarily, we 
find that besides oscillations of a' there are also oscilla- 
tions of the (/-shift with similar period and magnitude. 
Moreover, the periods of the oscillations show a depen- 
dence on the exchange splitting that is different from that 
of Ref. 12. 



II. SCATTERING BY A SQUARE WELL 

We consider a trilayer, N/F/N, consisting of a ferro- 
magnetic layer of thickness D imbedded between nor- 
mal metal layers, each of infinite thickness. The conduc- 
tion electrons move in a spin-dependent potential U a (x) 
where x is the coordinate perpendicular to the layers, 
and a specifies the orientation of the spin relative to the 
direction of the magnetization of the ferromagnet. Wc 
assume that the electrons in the F layer are described 
by a single-band Hubbard model 15 . In the Hartree-Fock 
approximation, this model yields one-electron energies 
£k,a = £k + I < n- a > where I is the on-site screened 
Coulomb interaction. To simplify the calculations, wc 
shift the energies in the F- region by —I < n |>. This 
brings the zero of energy to the bottom of the minor- 
ity band. In this case, the spin-up conduction elec- 
tron sees an attractive square well potential of depth 
A = /(< n t> — < n j>), a quantity defining the ex- 
change splitting between the majority and minority spin 
bands. The down- spin electron moves in a potential that 
is constant and equal to zero for all x. 

Thus, Ui(x) = for -00 < x < 00, Uf(x) = -A 
for -D/2 < x < D/2,Uf(x) = for x < -D/2 and 
x > D/2. 

If the zero of energy is taken at the bottom of the ma- 
jority band, the up-spin electron sees a potential barrier 
of height A, while the down-spin electron moves in con- 
stant zero-value potential. This choice of potential has 
been previously used by Bruno 10 ' 16 . Mills 12 assumes that 
the potential energy of the up-spin electron is Vq above 
the zero of energy (set at the bottom of the conduction 
band) whereas the down-spin has potential energy Vo+A. 
Thus, Bruno's potential follows by setting Vb = and our 
choice corresponds to taking Vq = —A. The fact that the 
results 12 for the D-dependencc of a' arc not very sensi- 
tive to the choice of Vq provides some justification for the 
present choice of the potential. Our main motivation for 
this choice is that the ^-component of the electron wave 
vector for spin up stays real (see Eq.(9)). This simplifies 
the summation over the transverse modes in Eqs. (4) 
and (5). 

The single-band Hubbard model implies that, owing to 
the high density of states, the effective mass of the itin- 
erant electrons in the ferromagnetic layer is larger than 
the free electron mass. This complicates the scattering 
problem since, besides the potential well, there is also a 
"mass barrier" in the F- region. Similar to Refs. 10 and 
16, we replace this position dependent mass by a con- 
stant average mass. It should be pointed out that an 
enhancement of the effective mass can also take place in 
normal layers due to electron-electron interactions. This 
effect is especially pronounced for Pd and Pt layers 13 . 

We first consider the reflection amplitudes rt,L Since 

the transverse momentum is conserved, we have rt,l = 

fm^m,n where = i"o*\kj^) where is the com- 
ponent of the wave vector along axis x. For our choice of 
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the potential, we have ^(fcjj = 0' and 16 



1 - exp(2ifc' D) 



1- [rL(fcl)] 2 exp(2ifc^) 



where 
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k± + k]_ 
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where k± and k\_ are the wave vectors in the N and 
F-regions, respectively. Since the expressions (4) and (5) 
are to be evaluated at the Fermi energy, we have 



and 
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where ku is the magnitude of the in-plane wave vector 
fen, and S 2 = 2m A/ft 2 . Eq. (9) implies that k]_ stays 
real over the whole range, (0, kp), of the variable ku. 
Next let us consider the transmission amplitudes. The 

reflection symmetry of the potential implies t m a n = tm]n- 
Due to the conservation of transverse momentum, we 
have tm]n = to"\k^)6 m>n where ^(^1) = 1 an<1 
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III. PARAMETERS ^,,4, 

We are now ready to evaluate the quantities A r and 
Ai. The transverse-mode sums in Eqs. (4) and (5) can 
be written as a sum over the in-plane wave vectors fe||. 
Converting the sum to a two-dimensional integral and 
using Eq. (9), we have 
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where L is the lateral dimension of the film. Owing to the 
presence of the term exp(2ik]_D) in Eqs. (6) and (10), it 
is convenient to integrate over the variable k ± . According 
to Eq. (9), the integration limits for this variable are 



fci = (k F + 5 2 )i 
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Using this relation and the fact that r| n „ = 0,4„ = l, 
equation (4) can be simplified to the following form 

A r = ]T(1 - Re4»,J (15) 

m.n 

Performing the summation with use of the prescription 
(11) and using Eq. (10), we obtain from Eq. (15) 

Ar = 2^ J k2 * A 
J(^) 2 + (rL) 4 ][l-cos(2fcl^)] 



l-2(rL) 2 cos(2fcl^) + (rL) 4 



(16) 



where 
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Substituting (r^ n = 0,t^ n = 1) into Eq. (5), we 
obtain 
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This result can be expressed using Eqs. (10) and (11) 
as an integral over k\_ 

r2 cfei 
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Making a substitution u = 2Dk\ , we obtain from Eqs. 
(16) and (17) 
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where 



r(u 1 ,U 2 ) = 



duu 



(1 - cosw){w 4 + [{u 2 - «|)3 - u} 4 } 
2wf(l — cosu) + l&u\u 2 (u 2 - u 2 ) 
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According to Eqs. (12) and (13), the dimensionless 
parameters u\ and u 2 are given by 



ui = 2Dh = 2D(k 2 F +5 2 )z 



(22) 



k 2 = S (13) 

Using the unitarity property of the scattering matrix, 
we have 



V m,nJ 1 V m,n/ 



l 



(14) 



u 2 = 2Dk 2 = 2DS (23) 
With the same substitution, Eq. (19) yields 
L 2 S 2 

A i = -—I i (u 1 ,u 2 ) (24) 
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where 
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Integrals (21) and (25) are evaluated numerically as a 
function of the parameter u 2 = 2DS by setting U\/u 2 = 
2.38. This value results from Eqs. (22) and (23) by taking 
Ef = 7cV and A = 1.5eV as adopted by Bruno 10 for the 
Co/Cu/Co system. We note that this set of parameters 
corresponds to replacing the x-dcpcndcnt electron mass 
by the free-electron mass. 

The results are plotted in Figs. 1 and 2. The data are 
given for the range < 2D5 < 60. Taking Uf — ir/a 
where a is the lattice constant of the ferromagnetic film, 
the upper limit of this range corresponds to D max ~ 
20a. From Fig. 1 we see that, for large values of D, 
the integral I r (D) tends to a constant equal to 0.08. As 
I r {D) approaches this value, it exhibits large quantum 
oscillations with an initial amplitude of 0.12. The period 
of these oscillations is about tt/S. This result can be 
traced to the singular behavior of the integrands of Eqs. 
(21) and (25). In fact, a numerical plot of the coefficient 
of the factor 1 — cosu in the integrand of Eq. (21) shows 
a sharp peak at u 2 — 2DS. Also the coefficient of the 
factor sinu in the integrand of Eq. (25) shows a peak 
at this value of u 2 . This behavior of the integrands is 
also responsible for the fact that the integrals (21) and 
(25) are insensitive to changes of the ratio u\ju 2 . For 
instance, on increasing this ratio from 2.38 to 10, the 
plots of the I r and Ii functions shown in Figs. 1 and 2 
sustain less than a 10% change. 

The amplitude envelope of the oscillations of the func- 
tion I r decreases slowly with increasing D as a power law 
D~ n where the exponent n ~ 0.6 ± 0.1. 

The integral (25), displayed in Fig. 2 as a function of 
the parameter 2d5, exhibits oscillations of similar char- 
acter taking place about Ii{D) = 0. The position and 
magnitude of the first peak of the function Ii (D) roughly 
agree with those of the oscillatory part of the I r (D) curve. 



IV. GILBERT DAMPING 

The excess damping constant a' is given by the second 
term on the right hands side of Eq. (3). Assuming that 
the increment the gyromagnctic ratio is small, we set 
7/70 ~ 1 and obtain with use of Eq. (20) 
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where M s is the saturation magnetization of the ferro- 
magnetic layer. 

Invoking Fig. 1, we see from Eq. (26) that, for 
D 3> Oi'(D) involves a D^ 1 dependence similar to 



previous works 1,4 ' 13 . Superimposed on this dependence 
are large slowly decaying short-wave oscillations with a 
period close to n/5. In the limit of D 1/5, the N/F 
and F/N interfaces are decoupled, and presumably con- 
tribute to a' independently 4 ' 13 . In Ref. 13, the interfaces 
were represented by monolayers acting as delta function 
scatterers. The expression for a' obtained in this ap- 
proximation is, to leading order in J sc i, proportional to 
J 2 d oc A 2 . This contrasts with Eq. (26) showing a linear 
dependence on A. A comparison of this result with the 
dynamic RKKY theory of Mills 12 is made in Sec. VI. 

Let us now estimate the magnitude of a' for the 
Cu/Co/Cu system. Consistent with the choice u\/u 2 = 
2.38 made in the evaluation of the integrals (21) and 

(25) , we take S 2 ~ 2mA/h 2 where A=1.5 eV. Letting 
M s ~ 1400 G, D ~ 40A and I r (D) w 0.1, we obtain 
from Eq. (26) a' ~ 3 x 10~ 4 . This is possibly an un- 
derestimate since a free-electron mass is assumed 10 . An 
enhancement of the average mass due to high density of 
states in the Co layer can increase the value of S 2 in Eq. 

(26) yielding a' of order 10~ 3 . 

Heinrich ct al. 17 studied, using FMR, epitaxially stabi- 
lized ultrathin fee films of Co deposited on bulk Cu(001). 
For a 10 monolayer (ML) Co film, they find a « 10~ 2 . 
Schreiber et al. 18 measured Gilbert damping of single 
crystal Fe x Co\- x alloy films prepared by sputtering on 
MgO(OOl) substrates. For pure Co films 200A thick, they 
measure a w 10 -2 for the hard and a~6x 10~ 3 for the 
easy direction, respectively. On the other hand, the pure 
Fe films show substantially lower damping (a « 2 x 10~ 3 ) 
that is practically isotropic 18 . On comparing our pre- 
dicted damping with the experimental results on Co 
films 17,18 , it is evident that the magnitude of the oscilla- 
tory part of a is within the error margins of the measured 
damping. 

Another way to estimate the Gilbert damping constant 
is to explore the current-induced switching of magnetic 
moments in Co/Cu/Co systems. Katine et al. 19 used 
thin-film pillars containing a thin (25A) Co film cou- 
pled to a thicker (100A) Co film. During the rotation 
of the magnetization of the thinner film, the magnetiza- 
tion of the thicker film remains fixed. This enables the 
determination of the polarity of the current bias asso- 
ciated with the spin-transfer excitations in the thinner 
layer. Moreover, the thicker layer acts as a spin sink. 
In this way, problems with spin accumulation in the Cu 
layers can be circumvented . The Gilbert damping con- 
stant is determined starting from the Landau-Lifshitz 
equation augmented by the current-induced spin transfer 
torque 20 . The stability of the given magnetic configura- 
tion of the pillar is determined by a competition between 
the spin-transfer torque and the Gilbert damping torque. 
From the way the critical switching current scales with 
the effective magnetic field, Katine et al. 19 find that an 
agreement with their experiment is obtained by assuming 
a fa 7 x 10~ 3 . This value is in good agreement with the 
FMR data 17 ' 18 . The current-induced switching has been 
also studied by injecting the current into the multilayer 
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through a point contact 21 . Interestingly, the damping 
constant deduced by this method is 10-50 times larger 
than the value obtained in the pillar geometry 19 . This 
could be related to the fact that, in the point-contact 
study, the spin-wave excitation is induced in a localized 
region exchange coupled to the unbounded ferromagnetic 
film. 

Due to the relatively small FMR linewidth of Fe 
films 18 , multilayers involving Fe layers seem more suit- 
able for the observation of the quantum oscillations. Ur- 
ban et al. 2 studied FMR on two Fe layers separated by 
nonmagnetic Au spacer. The thickness of the thinner 
layer is varied from 8 to 31ML. The thicker, 40ML, layer 
exhibits precession of negligible amplitude and serves as 
a spin sink for the spin currents emitted form the pro- 
cessing thin layer. The FMR linewidth follows a D^ 1 de- 
pendence without superimposed oscillatory component. 
For the 16ML film, the measured excess Gilbert damp- 
ing constant, a' « 2x 10~ 3 , is comparable to the intrinsic 
damping in the single Fe film 18 . 

Using Eq. (26), we estimate the theoretical value 
a' th for the 16ML Fe film imbedded between Au layers. 
We assume that the electronic parameters are similar to 
those of the Fe/Cu/ Fe system discussed by Hood and 
Falicov 22 . The conduction electron mass is taken, by 
these authors, independent of the material and the spin 
orientation and equal to m* w 4 x free electron mass. 
The exchange splitting for Fe is 22 A « 2.5 eV. The 
resulting value of the parameter S 2 — 2m*A/h 2 is 
2.6 x 10 16 cm~ 2 . Consequently, the period of the oscil- 
lations of the functions I r and Ii is n/5 w 2 x 10~ 8 cm. 
Since this period is comparable to the lattice constant 
of Fe, films differing in thickness by 1ML should be in- 
vestigated to detect the oscillations. Using this value of 
<5 2 , and taking M s w 1700G, D « 46l and I r « 0.1, 
we obtain from Eq. (26) the excess damping constant 
a' th « 1.7 x 10~ 3 . This is comparable with the experi- 
mental value 2 . 



From Figs. 1 and 2 we see that the oscillatory part of 
I r is of similar form and magnitude as Ii. This implies, 
that an observation of quantum oscillations of a' should 
be accompanied by an oscillation, of similar magnitude 
and period, of the ratio Ag/g . 

Fig. 2 implies that for thicker films, such that quan- 
tum coherence is suppressed, the shift of the (/-factor 
should be negligible compared to the measured value of 
a'. This prediction is in disagreement with the FMR 
data on Pd/Py/Pd and Pt/Py/Pt where Py is a permal- 
loy film 1 . These data yield a damping parameter a' that 
is proportional to 1/D and does not exhibit any superim- 
posed quantum oscillations. This suggests that the tri- 
layers are in a quantum-incoherent regime for which the 
present theory predicts Ag/ g^ ~ 0. In contrast, the ex- 
perimental value of this quantity is comparable with the 
experimental damping constant a' . For instance, for the 
Pd/Py/Pd system with D = 50l, a' « Ag/g w 10~ 2 . 
Thus we have 



Ag_ 

9oa' 



(29) 



It is interesting that an order of magnitude agreement 
with this result can be obtained by evaluating the ratio 
Ag/ (goce') from the model of two monolayers represented 
by a delta function potential 13 ' 14 . For Pd and Pt, the 
Coulomb interactions between the electrons need to be 
taken into account. In Ref. 14 we have shown that the 
effect of these interactions is to enhance the delta func- 
tion potential of the monolayer by the Stoner factor Se- 
As a result, the quantity Ai for the monolayer is changed 



to 
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V. GYROMAGNETIC RATIO 

Next we consider the renormalization of the gyromag- 
netic ratio 7. Putting 7 = 70 + A7 where A7 <C 70, we 
obtain from Eq. (2) 
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where Ag = g — go is the change of the g-factor due to 
the spin-pumping mechanism. 

Using Eqs. (20), (24), (26) and (27), we evaluate the 
ratio 



where /3 « kpJsdl '(2ef)- In what follows, we assume 
that J sc i is neg ative 23 . The origin of this antifcrromag- 
netic s-d interaction is quantum-mechanical mixing be- 
tween conduction and localized <i-orbitals. In Sec. II we 
found that the single-band Hubbard model 15 yields a fer- 
romagnetic interaction between the itinerant electron and 
the magnetization. This results from the Hartree-fock 
approximation which linearizes the on-site Coulomb in- 
teraction of two electrons with opposite spin. Thus the 
energy of the up-spin electron is lower relative to the 
down-spin one. There arc no consequences of this di- 
chotomy for the quantity a'. However, the sign of Ag 
changes as the interaction goes from the antiferromag- 
netic to the ferromagnetic one. 

The quantity A r for the monolayer becomes 14 
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Setting J sd 
yield 



-£ F /10 and S E « 5, Eqs. (30) and (31) 
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Applying this single-monolayer result to the model of 
two monolayers located at the N/F and F/N interfaces 13 , 
the same result is obtained for the ratio Ag/(g a'). We 
see that the theoretical result (32) is in order of magni- 
tude agreement with Eq. (29). 



VI. COMPARISON TO DYNAMIC RKKY 
THEORY 

We now compare the results of Sec. IV and V to the dy- 
namic RKKY theory of Mills 12 . Let us begin by consid- 
ering the excess Gilbert damping given by Eq. (26). We 
note that the formalism of the spin-pumping approach is 
very different from the dynamic RKKY theory 13,14 . Thus 
we limit ourselves to comparing only the final results. 

Using the relation G' = ^Mga', we obtain from Eq.(26) 
the excess Gilbert damping constant 



G 1 
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We wish to compare this equation with the excess 
damping constant obtained using the dynamic RKKY 
method (see Eq.(29) of Ref. 12) 



AG = 



fi 2 B kHe F )kt(e F )f(D) 
8ir 2 hD 
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where f(D) is a dimensionless function which, for 
k^(ep)D 3> 1, approaches a constant not far from unity. 
This function exhibits quantum oscillations that con- 
sist of short-wave oscillations with period 2/[fcT(ei?) + 
fcJ-(eir)] superposed on long- wave oscillations with period 
2/[fct(ei?) — fc^(ei?)]. In contrast, the function I r (D) ex- 
hibits only short-wave oscillations with period tt/5. It 
approaches, as D increases, a constant close to 0.1. 

As far as the magnitude of the damping is concerned, 
the main difference between Eqs. (33) and (34) is that 
the factor (e f)^ (e f) is replaced by S 2 . This implies 
that the magnitude of G", for large D, is reduced com- 
pared to AG by a factor of order A/sf- A more serious 
discrepancy is that Eq. (34) yields a nonzero damping 
even when the splitting of the up and down-spin bands 
vanishes ( ^(sf) = ^(ef))- According to Eq. (4), A r 
vanishes in this case since the reflection and transmission 
amplitudes for the up-spin cancel those for the down-spin 
electron. This feature of the spin-pumping theory has 
been emphasized by Tserkovnyak et al. 4 when compar- 
ing their formula for excess damping with the result of 
Berger 3 . The requirement of nonzero spin-splitting also 



follows from the general dynamic RKKY approach 13 to 
Gilbert damping. We return to a possible resolution of 
this difficulty at the end of this section. 

In what follows, we argue that the absence of long- wave 
oscillations in I r (D) is an upshot of the simplified poten- 
tial (Ui(x) = 0). Similar argument has been raised by 
Bruno 10 for the interlayer exchange coupling. When both 
the up and down-spin electrons are scattered, a long-wave 
modulation of the short-wave oscillation is expected for 
the D-dependence of this coupling. 

To substantiate this idea, we derive A r for the case 
when both f7j.(x) and U^(x) are nonzero. Using Eqs. (6- 
10), we obtain from Eq. (4) 



T 2 f kp ( 
A r ~— dkfyiB^B^lB^B 1 (35) 

-[l-irini-irL) 2 }-^ 

x[l + (rl) 2 + (rU 2 + (rlrU 2 ] 
+{(rini-(rimi-(rl) 2 } 

x cos(2fc}£>) + {r^ <=► r^} cos(2fc[L>) 
-{rlrUirlf + (d) 2 ]} cos[2(fcl + k{)D] 
-{rlrUl + (rlrL) 2 } + (rl?i) 2 [l - (rjj 2 ] 
x[l-(ri) 2 ]}cos[2(/cl 

where 

= 1 - 2(r^) 2 cos(2k { ^D) + (r^) 4 (36) 
Generalizing Eq. (9) to both spin components, we set 



fc[ a) = [(4 a) ) 2 -# 



(37) 



where k^ = (k'p + <5 2 )2 is to be identified with 
k^(e F ) of Ref. 12. 

To determine the periods of the oscillations of A r (D), 
we consider the asymptotic estimate of the integral (35) 
for D large. We note that the major contribution to 
this integral arises from the immediate vicinity of the 
end points of the integration interval (0, kp) and from 
the vicinity of the stationary points of the function k±"\ 
From Eq. (37) we see that the derivative of this func- 
tion vanishes at the stationary point fc|| = 0. However, 
the contribution of this stationary point to the integral 
(35) vanishes owing to the factor fcii which causes the 
coefficients of all cosine terms to vanish. 

For the same reason, the partial integration method 
yields a nonzero contribution only from the end point 
fcii = Uf- This conclusion holds also for the simplified po- 
tential used in Sec. III. In fact, setting = 0, Eq.(35) 
reduces to the integral (16) which is dominated by vicin- 
ity of k\{k\\ = kp) = S. 
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Consequently, Eq. (35) yields three short-wave oscil- 
latory terms with periods 7r/<5f, tt/S±, 7r/(5| + Si) and a 
long-wave term with period 7r/(<5j — Si). Hence, once 
electrons of both spin orientations are allowed to scatter, 
the spin-pumping approach exhibits oscillations similar 
to those of Ref. 12. The main distinction is that the 
wave vectors k^(ep) determining the periods of the os- 
cillations in the dynamic RKKY theory 12 are replaced by 
the wave vectors S a . 

We now return to the aforementioned problem con- 
cerning the product fc T (e F )fc^(£ F ) in Eq. (34). First, we 
examine the integral (35) for the case when S^ = Si (im- 
plying = r^). In this case, it is possible to verify 
explicitly that A r {D) — > as D — » oo. Specifically, the 
rapid oscillations tend to cancel the contribution of the 
three short-wave cosine terms in the integrand. The re- 
maining terms, upon substituting = r^, exactly can- 
cel so that the integral (35) vanishes. This is in accord 
with the more general formula (4). Moreover, it provides 
an independent check that the rather complicated inte- 
grand of Eq. (35) is correct. Based on this result and the 
fact that, upon letting = 0, this equation goes over 
to Eq. (16), we propose the following generalization of 
Eq. (20) 

A r (D) oc L 2 ((5 T - Si) 2 F r (D) (38) 

where F r (D) tends to a constant for D large and con- 
tains a superposition of slowly decaying oscillations with 
periods 7r/(5|, 7t/<5|, 7t/(<5| + <5|) and 7t/(<5| — <5j). Eq. (38) 
also suggests that replacing the product k' (e p)k^ (e f) 
by [k^(ep) — fc|(ei?)] 2 may endow Eq. (34) with a cor- 
rect dependence on the splitting of the up and down-spin 
electron bands. 

In the preceding section we came to the conclusion 
that quantum oscillations of the damping constant a 1 
are accompanied by oscillations of similar magnitude and 
period of the relative shift, Ag/g , of the gyromagnetic 
ratio g. This spin-pumping result disagrees with that 
of the dynamic RKKY method 12 . In this method, the 
rcnormalization of g stems from the frequency deriva- 
tive, called Ai, of the transverse dynamical susceptibil- 
ity. Though Ref. 12 points out the general importance 
of this quantity in the spin dynamics, it reports rough 
estimates indicating that the renormalization effects due 
to Ai are small. 



VII. DISCUSSION 

In the itinerant model of a bulk fcrromagnet, the g- 
shift and the damping a' have been attributed to the 
combined effect of spin-orbit interaction and the electron- 
lattice interaction 24 . Mizukami et al. 1 assume that Ag 
and a' observed on Pd/Py/Pd and Pt/Py/Pt trilayers 
are due to similar interactions taking place at the inter- 
faces. 



The results outlined in Eqs. (30)-(32) represent a 
strong departure from these ideas. In the framework 
of the spin-pumping theory, the enhancement of the gy- 
romagnetic ratio comes from the angular momentum of 
the conduction electron spin carried away from the F- 
film into the N-layers. It is the spin current that is in 
the phase with the precession that contributes to the g- 
factor. The out of the phase part gives rise to the excess 
damping a'. 

Now in order that this mechanism works, it is neces- 
sary that the conduction electron spins are at equilibrium 
with the normal metal reservoirs. Otherwise, spin accu- 
mulation takes place which blocks the flow of the spin 
current 4 . The equilibrium is established by spin-lattice 
relaxation in the normal metal 25 . This relaxation pro- 
cess involves the spin-orbit interaction combined with the 
electron-lattice interaction similar to the mechanism of 
Ref. 24. In the spin-pumping approach, the conduction 
electron spin is first transported into the normal metal 
where it relaxes. Thus we may view the relaxation of 
the ferromagnetic magnetization, in this approach, as a 
nonlocal process mediated via spin current. On the other 
hand, the itinerant electrons of the bulk ferromagnet re- 
lax via local spin-orbit interaction 24 . Clearly, the spin- 
orbit interaction is an inevitable component for magne- 
tization relaxation in both mechanisms 4,24 . 

It remains to address the reasons for the absence of 
quantum oscillations in the FMR experiments 1 ' 2 which 
study the Gilbert damping and the (/-factor as a function 
of the thickness of the ferromagnetic film. It is, of course, 
possible that another mechanism [see Ref. 1] intervenes 
and obliterates the predictions of Eqs. 26 and 27. 

Another possibility is that the oscillations are washed 
out by interface roughness. Note that the wavelength 
of the short-wave oscillations, considered here, is about 
b°A for Cu/Co/Cu and only about 2°A for Cu/Fe/Cu. It 
is conceivable that roughness introduces fluctuations of 
D larger than these wavelength. Averaging over these 
fluctuations strongly suppresses the amplitude of the os- 
cillations. Consequently, one essentially sees an average 
which is a constant of order 0.1 for Fig.l, and zero for 
Fig. 2. This is in apparent disagreement with the mono- 
layer model invoked in Sec. V to explain the large shift of 
the (/-factor observed on Pd/Py/Pd systems 1 . To clarify 
this point, we note that the monolayer model assumes 
that the quantum oscillations in the ferromagnetic film 
rapidly decay with the distance from the interface. In the 
extreme case of overdamping, the function h{D) does not 
change sign upon increasing D by one wavelength form 
the origin. Thus fluctuations of the thickness do not av- 
erage the g-shift to zero in the monolayer model. 

Though short-wave oscillations are expected to be 
washed out by surface roughness, there are also long- 
wave terms in the expressions for the Gilbert damping 
and the shift of the gyromagnetic ratio that may survive 
the averaging over D . As seen from Eq. (35), these terms 
appear in A r when electrons of both spin orientations are 
subject to scattering. It can be shown in a similar way 
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that these terms are also present in the expression for 
Ai. Since the period of these oscillations is n/(5^ — <5jJ, 
multilayers with small exchange splitting should be used 
to increase the chance of observing the long-wave oscil- 
lations. 

An alternative way of observing the quantum oscilla- 
tions is suggested by the fact that I r and Ii are functions 
of the product 2DS. Thus, instead of changing the thick- 
ness of the ferromagnetic film, we may vary the param- 
eter 8 while keeping D fixed. For the case when both 
[7f and are nonzero, we have <5 2 = 2mA a /h 2 where 
A a is the depth of the potential well for electron of spin 
a. These potential wells can be varied by shifting the 
chemical potentials, fi a , of the electrons in the N-region. 

One possibility is to apply a potential difference across 
the N/F interface. This would change the parameters <5j 
and 6 1 by the same amount. Consequently, only short- 
wave oscillations could be tuned in this manner. It is 
worth noting that, for the N/F interfaces considered here, 
a voltage drop of order leV is needed to sweep through 
several periods of the short-wave oscillations. To mini- 
mize Joule heating associated with such a voltage drop, 
the interface resistance should be large. This could be 
achieved by a thin layer of oxide incorporated into the 
interface. 

Tuning of long-wave oscillations requires that the dif- 
ference A | — A i is varied. It is tempting to consider the 
shift A/i = (pi — Hi) due to the spin polarized current 
normal to the interface 26 . Using Eq. 12 of Ref. 26, we 
estimate that the current density required to generate a 
shift of leV is of order 10 10 A cm~ 2 which rules out this 
idea. 

It is of interest, both experimentally and theoretically, 
to examine other processes contributing to the decay of 
multiple reflections in ideal multilayers. Recent work 
of Stiles and Zangwill 27 is of relevance in this context. 
These authors show that the transverse [i.e. perpen- 
dicular to m] component of the spin current that flows 
from the normal metal into the ferromagnet is absorbed 
within a distance of several lattice constants. In the spin- 
pumping approach 4 , it is also the transverse spin current 
that determines the Gilbert damping and the gyromag- 
nctic ratio. The parameters A r and Ai give the mag- 
nitude of these currents. Now, the present calculation 



of these parameters shows slowly-decaying oscillations in 
disagreement with Ref. 27. Also the numerical results of 
Mills 12 do not show fast decaying oscillations of Gilbert 
damping. We note that Stiles and Zangwill 27 have iden- 
tified three distinct processes contributing to the absorp- 
tion of the transverse spin current: (1) spin dependent 
reflection and transmission; (2) rotation of reflected and 
transmitted spins; and (3) spatial precession of spins in 
the ferromagnet. Out of these processes, it is the first 
one that the present theory and that of Ref. 12 takes 
into account. It should be interesting to incorporate the 
other two processes into these theories to see if a fast 
decay of the quantum oscillations ensues. 

It should be emphasized that these decay processes 
do not affect the longitudinal component of the spin 
current . Incidentally, it is also this component which 
determines the exchange coupling between ferromagnetic 
films separated by normal metal spacer. Thus, we ex- 
pect that the -D-dependence of the the coupling exhibits 
slowly decaying oscillations. This agrees with the theory 
of Bruno 10 and its experimental verification 11 . 

In view of the predicted fast decay of the transverse 
spin current, the independent oscillations may be decay- 
ing faster than calculated in the present simplified model, 
and therefore less easy to detect than the oscillations of 
exchange coupling. In this context, it is worth noting 
that the decay of the transverse spin current depends on 
the Fermi surface mismatch between the ferromagnetic 
and normal metal Fermi surfaces 27 . A slower decay is 
predicted when the mismatch is small (k F /kp w k^ /kp). 
From this point of view, multilayers with small exchange 
splitting are preferable in order to mitigate the effect 
of the decay processes 27 on quantum oscillations of the 
Gilbert damping and gyromagnetic ratio. 
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FIG. 1: Integral I r (ui,U2) defined in Eq. (21). The plot is 
obtained by numerical integration as a function of U2 = 2SD 
for u\/u2 = 2.38. D is the thickness of the ferromagnetic 
layer, and S = (2mA/S 2 )^ 2 where A is the depth of the 
potential well. 



FIG. 2: Plot of the quantity /i(«i, 1*2) evaluated from Eq. 
(25) as a function of 25D for u\jui = 2.38. 
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